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Abstract 

We calculate the lepton electric dipole moments in a class of supersymmetric seesaw models and 
explore the possibility that they may provide a way to probe some of the CP violating phases 
responsible for the origin of matter via leptogenesis. We show that in models where the right 
handed neutrino masses, M/j arise from the breaking of local B-L by a Higgs field with B-L=2, 
some of the leptogenesis phases can lead to enhancement of the lepton dipole moments compared 
to the prediction of models where M/j is either directly put in by hand or is a consequence of a 
higher dimensional operator. 
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I. INTRODUCTION 



A presently popular way to understand the origin of matter-anti-matter asymmetry in the 
universe is to start with a mechanism to generate lepton asymmetry in the early Universe 
using the CP violating effects in the decay of a heavy right handed Majorana neutrinosj^ 
and let the sphaleron interactions above the electroweak phase transtion temperature 
convert them to a baryon asymmetry. This is possible since the sphaleron interactions which 
violate B+L symmetry are in thermal equilibrium above the eletcroweak phase transition 
temperature. The right handed neutrino decay is of course not the only way to generate the 
pre-electroweak lepton asymmetry. Several other ways have been discovered for the purpose 



right handed neutrinos since they also provide a simple way to understand small neutrino 
masses via the seesaw mechanism jsl. In particular, the existing experimental information 
on neutrino oscillations via the seesaw mechanism effects the nature of the right handed 
neutrino mass pattern, which in turn effects the magnitude of the baryon asymmetry. It is 
quite interesting that detailed analyses that use the exisitng neutrino data do indeed give 
the right magnitude for baryon asymmetry as well as important insight into the pattern of 
right handed neutrino masses 0]. 

As noted in the original pioneering work of Sakharov, CP violation is an essential in- 
gredient in the generation of any particle-anti-particle asymmetry. In the present case, the 
CP violating decays of the right handed neutrinos arise from the phases in right handed 
neutrino couplings and we will call them leptogenesis phases. It is clearly important to 
seek low energy manifestations of the leptogenesis phase for several reasons: first, this will 
improve our understanding of the right handed neutrino mass matrix which plays a crucial 
role in the neutrino mass physics at low energies; secondly, it may shed light on the origin 
of the seesaw mechanism, which will then provide a useful window into physics beyond the 
standard model. Moreover, since the seesaw formula provides some connection between the 
low and high scale phases in the theory, understanding leptogenesis phases may be a guide 
to the CP violating effects in neutrino oscillations, with its many experimental ramifications. 

There have been a great deal of discussion of this issue in literature 0]. Specific models 
where the neutrino mass phase and leptogenesis phases are directly related have also been 
discussed in several papersjsl providing one way to probe the latter. However in all previous 




We will however focus in this paper on the decay of heavy 
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discussions of this issue, the masses of the right handed neutrinos Mji are either put in by 
hand or are assumed to arise from nonrenormahzable operators. It is however well known 
that in models that contain Higgs bosons A with B-L=2, arises from a renormalizable 
coupling /z/^z/^A. Examples of such theories are left-right symmetric models with triplet 
Higgs fields or SO(IO) models with 126 dimensional Higgs fields, SU{5)l x SU{5)ji models 
with Higgs fields in (1, 15) © (15, 1) representation. In such theories, there are new renormal- 
ization group running effects between the GUT (or Planck) scale and the seesaw scale. We 
showed in a previous paper that the presence of these renormalizable couplings can lead 
to quantitatively different effects in seesaw induced lepton flavor violation. In particular, 
the ratio B{fi e + 7)/i?(r — > /i + 7) depends on the SUSY breaking parameters in very 
different ways in different theories. 

In this paper, we calculate the electric dipole moments (edm) in the class of supersymmet- 
ric seesaw models with the / coupling and discuss how it helps in probing the leptogenesis 
phase. 

Since we are working within supersymmetric seesaw models, (i) we need to make assump- 
tions about the nature of supersymmetry breaking and (ii) the embedding of MSSM into 
new physics at the seesaw scale. 

As for the first point, in common with many discussions in the literature, we will assume 
that the TeV scale theory is the minimal supersymmetric standard model, which solves 
the gauge hierarchy problem plus the additional feature that the neutrinos have Majorana 
masses arising from the seesaw mechanism. As far as SUSY breaking is concerned, we will 
work within the context of minimal MSUGRA models and assume that around the GUT 
or the Planck scale, all soft SUSY breaking scalar masses are universal as are the gaugino 
masses and the A-terms are proportional to the Yukawa couplings. Furthermore we will 
assume that there are no overall phase in the A terms, which is guaranteed if for example, 
MSSM is part of a left-right symmetric SUSY model near the GUT scale^|. The resulting 
theory is a very economical one with only five parameters characterizing the complete susy 
breaking sector of MSSM and is therefore quite predictive. More importantly, the restrictive 
nature of this assumption means that at low energies, we are only tracking the effect of 
the phases in the right handed neutrino couplings responsible for leptogenesis and no other 
phase unrelated to neutrino physics is present. Once we give up these assumptions, new 
phases could come in and confuse the issue. 
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As far as the high theory theory is concerned, to implement our premise of having seesaw 
mechanism arise out of a B-L = 2 Higgs field, we will consider the high scale theory to be 
SU{2)i X [/(l)/^^ X U{1)b^l based with a standard model singlet A with B-L =2 breaking 
the B-L gauge symmetry and also giving the seesaw mechanism. This model can arise as 
an effective theory from left-right or SO(IO) theories. This extended model leaves the low 
energy predictions of MSSM uneffected. 

The weak scale values of the supersymmetry breaking parameters can then be derived 
by the renormalization group extrapolations using the standard techniques. Since it is the 
dimension four terms in the Lagrangian which are largely responsible for the running, how 
the right handed neutrino masses are generated in the Lagrangian does make difference to 
low energy phenomenology. We discuss the impact of this extra running effect on the lepton 
dipole moment observable which depends on the leptogenesis phases. 

The main result of our investigation is that if the right handed neutrino mass in the 
seesaw mechanism arises from the vev of a B-L=2 Higgs field, it has the effect of enhancing 
the lepton electric dipole moments over models where the right handed neutrino mass is put 
in by hand. We give two examples to illustrate this point: one where the neutrino mixings 
arise purely from Dirac Yukawa coupling of the right handed neutrinos and a second one 
where it arises from their Major ana coupling /. We also discuss the case of a 3 x 2 seesaw 
models. In all the B-L=2 Higgs field is responsible for righthanded neutrino masses. 

We present the results of our calculation for electron amd muon edms for the above 
models making sure that we stay in the range of parameters which fit all neutrino data and 
generate the right amount of lepton asymmetry. We find that leptogenesis phase can lead 
to enhanced edms for leptons in the presence of the / couplings than in their absence (i.e. 
where the RH neutrino has a bare mass term.), although the effects are still small. Most 
optimistic values for electron edms in our models are between 10^^^ to 10~^^ ecm and for 
muon edm between 10^^'' to 10"^^ ecm. It is encouraging that there are ideas and plans for 



drastic improvement in the search for the lepton edms in the future [1^ 3, 13- -^^^ instance, 
search for electrons edm upto the level of 10"'^^ ecmlldl and muon edm to 10~^^ ecm have 

n 

been contemplated |13l| . 

Although the values predicted by our analysis are small, high precision searches for lepton 
edms to the contemplated level can still teach us something useful. For instance, if the edms 
discovered are above our predictions, it will mean one of several things: either there is gross 
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departures from the assumed forms for supersymmetry breaking terms or if independent 
experiments such as at LHC or Tevatron have confirmed the SUSY breaking assumptions 
(universahty and proportionahty), we would have to conclude that baryogenesis does not 
originate via leptogenesis. Either of these conclusions would be very important in our search 
for physics beyond the standard model. 

We have organized this paper as follows: in sec. 2, we present a generic model based on 
the gauge group SU{2)l x [/(l)/^^ x f/(l)B-L to discuss the general calculational framework 
and the basic idea of our work; in sec 3, we discuss the renormalization group effects of 
the neutrino sector on the supersymmetry parameters, in sec. 4, we present a qualitative 
explanation of how the enhancement of edm can arise in our approach; in sec. 5, we briefly 
discuss the calculation of baryogenesis; in sec. 6 we discuss the first neutrino mass model, 
give the neutrino mass fits, calculation of lepton asymmetry and present the results for 
lepton edms; in sec. 7 and 8, we repeat the same discussion for two other models; in sec. 9, 
we summarize our results and present the conclusions. 

II. A GENERIC MODEL AND THE CALCULATIONAL FRAMEWORK 

To proceed with our discussion, let us write the superpotential for our model, which must 
be invariant under the gauge group SU{2)l x U{l)i.^j^ x U{1)b~l'- 



Here L, e*^, are leptonic superfields; Hu^, A, A are the Higgs fields. We do not display the 
quark part of the superpotential but it is same as in the MSSM. 

The first point to notice is that it is a very simple extension of the MSSM superpotential, 
that leads to a vacuum with <A> =< A >= w/j. This leads to the RH neutrino mass 
matrix Mr = fvR and clearly it arises from a renormalizable term in the Lagrangian. 

This superpotential has also another advantage that it solves the /i problemjl^ i.e. it 
predicts that yU ^ ^,3/2, where m^i2 is the superpartner mass, which is of order of a TeV. 
To see this, notice that this potential has an R-symmetry under which W — > —W\ (the 
fields transform as {Hu,e'^,d'^) — > i{Hu,e'^,d^)] {H^.v^.u'') —i{Hj^,v^,u^); S —S and 
all other fields remain unchanged under R-symmetry). The vanishing of the F-terms then 
imply that < 5* >= and < 5* > arises from supersymmetry breaking terms. As a result 
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/i ~< 5* >~ m3/2. 

The seesaw formula (type I) can be written as: 

= -Yjr^Y /^^^^'''^ (2) 

vr 

We can choose a basis where the couphng matrix / is diagonal as is the charged lepton 
Yukawa coupling matrix Yg. The neutrino Dirac coupling matrix Y,y is then a 3 x 3 complex 
matrix whose diagonal elements can be made real by suitable redefinition of the phases of 
the ii and fields. This leaves us with nine real parameters in the Yt, matrix and six phases. 
We rewrite the Y^, matrix as follows: 

Y^ = V*Y^W^ (3) 

where V = P1VP2 with = diag(e^'^^^ , e*'^'^'^ and V and W are CKM type matrices with 
only one phase in each. Thus all six phases are accounted for. 

To see which of the phases are responsible for lepton asymmetry, let us write down the 
formula for lepton asymmetry, n^jl^: 

n,^Y.MyX]hn^) (4) 

where symbol 1 in the above equation denotes the lightest right handed neutrino eigenstate. 
Note that using EqI21 one can conclude that only the phase in V and those in Pi contribute 
to lepton asymmetry. The phases in P2 completely drop out. Our goal will therefore be to 
see how we can measure the phases in Pi and V by low energy experiments. 

Secondly, we find two results for the three generation case that we state below. 



A. Independence of light neutrino mixings and leptogenesis phases 

The first point is that the lepton asymmetry parameter ne is independent of the PMNS 
mixing angles. To show this first note that 

= U*MtU^ = -Yjf-iY,mo (5) 

where mo = ^ and U is the PMNS mixing matrix. We can now invert the seesaw 

relation to conclude that 

Y, = tf/^0{Mty^^U^ (6) 
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where O is a complex matrix with the property that 00 ^ = 1[17\. The set of matrices 
O in fact form a group analogous to the complex extension of the Lorentz group. Using 
EqEl in Eq0J we see that lepton asymmetry is independent of the PMNS mixing angles 
in f/|il8.]. This result is important because what this means is that if there are three right 
handed neutrinos, low energy experiments such as neutrino oscillations and neutrinoless 
double beta decay give us no information about the leptogenesis phase. 



B. Vanishing of lepton asymmetry for degenerate neutrino masses 

To see this, we rewrite the formula for n£ in a slightly different way. Using the seesaw for- 
mula in Eg El and assuming that there is a mass hierarchy among the right handed neutrinos, 
we can write [l9| 

cx Im[Y,MlY^]n (7) 

Using this and EqlHl we can write 

ne oc Im[f/'0{Mtr0^f/' (8) 

From this we can conclude that when the neutrino masses are degenerate i.e. = ml, 
then ng = since fu is real. This result, to the best of our knowledge does not exist in 
literature. A similar result for the manifestly CP violating observable Pu^-ue — Pu^-ve 
noticed in Ref.j^], where it was pointed out that it vanishes for degenerate neutrinos. 

An implication of this result is that if neutrinoless double beta decay results imply m > 0.1 
eV(see e.g. then the lepton asymmetry from RH neutrino decay will be suppressed by 

an extra factor of 100 and we might then have to look for extra enhancements or completely 
different way of obtaining the baryon asymmetry of the universe. 



III. RGE'S AND LOW ENERGY EFFECTS OF LEPTOGENESIS PHASES 

Let us now proceed to the main part of our discussion, which involves the supersymmetry 
breaking superpartner proiperties at low energies. 

It is clear that one simple way to measure high scale parameters of a theory is to look 
for low energy observables that owe their origin or magnitude to the renormalization group 
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evolution (RGE) from the high scale to the low scale (i.e. without RGE effects, the low 
energy theory will predict the observable to be very tiny.) In our case such observables are 
lepton flavor violating transitions as well as lepton electric dipole moments (edms). In the 

4 

absence of the RGE effects, we expect ^ e + 7) to be proportional to ~ 10^^^ As 
far as the edm of leptons go, in the low energy theory without any other effects i.e. standard 
model with a Majorana mass for the neutrino, we would expect rfg — '^g^" to be of order 
10"^*^ ecm. 

In supersymmetric models, the presence of superpartners introduce new contributions to 
both LEV as well as edm observables at the one loop level, where virtual superpartners flow 
inside the loop. The resulting LEV and edm effects are therefore not only enhanced but they 
are also sensitive to the nature of SUSY breaking terms. In this paper we will assume a very 
restrictive form for the SUSY breaking terms which defines the so called mSUGRA models. 
According to this assumption, the scalar masses and all gaugino masses are universal at the 
GUT or the Planck scale; the A terms are proportional to the Yukawa couplings. These 
assumptions are motivated in several SUSY breaking scenarios; they not only reduce the 
number of parameters in the theory but may be necessary for understanding the observed 
suppression of flavor changing effects in both quark and lepton sectors. Eurthermore, we 
will assume that the overall phase in the A term is zero. This can happen for example, if 
the theory is part of SUSY left-right model at the GUT scale [l^]- This is necessary in order 
to resolve the so-called SUSY CP problem. 

The proportionality assumption coupled with the assumtion of no overall phase in the A 
term guarantees that all the CP violating phases in the lepton sector sector are in the Yi, 
and therefore any CP violating low energy effect after renormalization group evolution is 
only sensitive to the phases that are manifested in the leptogenesis as well as the neutrino 
oscillations via the seesaw mechanism. 

As mentioned, we need to extrapolate the SUSY breaking parameters from the GUT or 
Planck scale down to the seesaw scale where we will perform the calculation of the electric 
dipole moments of the electron and the muon. This extrapolation will always involve the 
couplings Yi, and the diagonal coupling matrix Y^. Extrapolation below the seesaw scale 
only involves the Yg which is mostly diagonal to leading order and does not effect any of our 
results. We will keep only one loop effects. The relevant renormalization group equations 
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are: 



fjv Y 1 

i = -^[{ml + 2mlm^Ye + {ml + 2Hl)YX + ^Y,^^lYe + Y,^YemllO) 



dm 



dt 167r2 



,.,2 1 



Let us also write down the renormalization group equations for the A parameters- specifically 
the A^. 

^ = jt^lMTriSYjY, + Y.^Ye) + W}Y, + Y^ (H) 
- -gl- 3/2^Ll] + Ye[Tr{6AaY^ + 2^,^/) 
+ 4y,U, + 2yJ-A, + 6g'^M2 + 2glMR + 3^^_^Mb-l]] 

2 [A.[Tr(3y;ty:, + rjy;) + 5yX + y^y^ (12) 



dA. 



dt IGtt 

+ -3^2' -gl- ^I'^gl-L + + n[rr(6^„Fj + 2^.f;) 

+ 4yJ-A, + 2Y}Ae + 6^2'M2 + 2glMR + 3^|_^Mb-l] + 8A//ty,] 

While in presenting actual results, we numerically integrate these RGEs to obtain their 
contribution to physical observables such as lepton edms, it is instructive to give some ap- 
proximate analytic expressions for the different results to get feeling for our final predictions. 
We do this in the following section. 



IV. ESTIMATES OF THE ELECTRON AND MUON EDMS AND CONNEC- 
TIONS TO LEPTOGENESIS PHASE 

To see the new effects from the / couplings in a qualitative manner, we note that the 
dominant supersymmetric contributions to the lepton dipole moments come from the one 
loop diagram involving the bino. A rough estimate of this is 

on ivi^^^ 

If there is an overall phase in A(_, then one can easily estimate d^ ~ 10~^^sin5 ecm. To be 
consistent with present limits i.e. < 1.7 x 10~^^ ecm, one must have phase 5 < 10~^. This 
is the so called SUSY phase problem. We assume that in the theory the SUSY CP problem 



has been solved, (e.g. by embedding the MSSM into a SUSY LR model as mentioned) so 
that we can set 5 = 0. In this limit, in the absence of the RGE effects, since is real, 
in this limit, all edms vanish. If this term dominated the contribution to the lepton dipole 
moments, then since the A term id proportional to charged lepton masses, we would have 
the scaling law d^^/dl = m^/mf.. In any theory where this law holds, the present upper limit 
on electron edm of < 4.3 x 10~^^ ecm would imply that (i^ < 8.6 x 10~^^ ecm. Also this 
scaling law is a signature of the simplicity of the underlying theory. In the supersymmetric 
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context, there exist exceptions to this scaling law 

After RGE effects are taken into account, lepton edms receive contributions from Y^Y^ 
which is not a diagonal matrix as well as /^/ and Y^ which are diagonal matrices at the GUT 
scale. The and Ai at the seesaw scale are polynomials in these matrices and in order 
to get a nonvanishing contribution to c?^, we must look for imaqginary diagonal elements in 
the product AfM|, which will involve products of Y^Yy, f^f and Y^. 

The leading order term in A^M| in the absence of the / couplings (i.e. a bare mass term 
for Ufj), which has a complex 11 entry is of the form 

d', (X Im[moYeYXYi\YXy]. (14) 

From Eq. (3), it is easy to see that YjY^ is independent of the matrix V that is responsible 
for leptogenesis. One might conclude from this that simple RGE extrapolations like this 
term will never enable us to probe the matrix V. However, as has been argued inj^l, since 
YJYi, is sensitive to the unitary matrix W in Eq. (3). So in principle, one can get W from 
the experimental information about supersymmetry breaking sector. Then, if one is given 
the right handed neutrino mass spectrum and the form of Y^, one can get all the matrix 
elements of V using the relation 

v^fv = -{Y;^y^w^u*Mtu^w{Y;})-^ (15) 

and hence the baryogenesis phases. This would require knowledge of all the matrix ele- 
ments of slepton mass matrix. While this is possible in principle, one would like a more 
"experimentally feasible" way to probe the phases in V . This is what we pursue here. 
There are two ways in which the RGE extrapolation can still leave traces of the matrix 



V . One way discussed in 



24| is the following. In the case where the right handed neutrino 



masses are hierarchical, the simple RGE expression in the above paragraph gets modified by 
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threshold effects as follows: the simple expression in Eg 1141 remains valid above the highest 
right handed neutrino mass, M3; once we go below this i.e. for M2 < /i < M3, only the 
second and the lightest RH neutrinos contribute. This means that in the matrix Y^,, only 
Yjj2i and Yy^u contribute and below M2 only Yy^n entries contribute. Therefore traces of 
the matrix V remain in the low energy expressions. To get a large effect, a significant gap 
between the M3 and Mi is required. 

Another way that V can be visible is the presence of the / terms, which give the leading 
contribution to edm to be Y^YJY^ p f'Y^ instead of what is given inHH Thus in this case, 
there is a higher order effect even in the absence of threshold effects. 

Note further that each time a term of the form Y^Y^ or f'^ f appears, there is a loop 
suppression factor eioop — Yq^^^ (1?^) ~ ^■^^ ~ ^^^^ Taking this into account we see that 
the second contribution (i.e. one in the presence of the / couplings) is enhance by the inverse 
of one loop factor compared to the first one. Also, there is an additional effect coming from 
the gap between the Mjj or Mpi and M3. Therefore barring unexpected cancellations, we 
expect the edm contribution in the presence of / contributions to be bigger. In our numerical 
integrations to present the final numbers for edms, we keep the threshold effects. 

V. OVERVIEW OF BARYOGENESIS ESTIMATE 

As already mentioned earlier, we use the mSUGRA framework to calculate the electric 
dipole moments of the muon and the electron. The value of the universal scalar mass mg, 
universal gaugino mass mi/2, universal trilinear term Aq, tan/5 and the sign of mu as free 
parameters determine our final result. We also assume that there is no phase associated 
with the SUSY breaking. The Yukawa and/or the Majorana couplings are responsible for 
CP violation in these models. 

The mSUGRA parameter space is constrained by the experimental lower limit on nih and 
the measurements of 6 ^ 57 and the recent results on dark matter relic densityj^]. For 
low tan/?, the parameter space has lower bound on mi/2 from the light Higgs mass bound 
of > 114 GeV. For larger tan/5 the lower bound on mi/2 is produced by the CLEO 
constraint on the BR(6 — > 57). The recent WMAP results have led to stronger constraints 
on dark matter density which in turn reduces the allowed parameter space of mSUGRA 
mostly to the co- annihilation region for mo, mi/2 < 1000 GeV. 
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We calculate the baryon to photon ratio ?7_b in three different models for neutrino masses. 
The SUSY breaking parameters do not play a role in this discussion. The lepton asymmetry 
is produced from the out of equilibrium decays of right handed neutrinos and this gets 
converted into baryon asymmetry through the sphaleron processes. The formula for rjs is: 

Vb = — = ag'i}_^— (16) 

g* 

where a relates the B-L asymmetry to baryon asymmetry and ^ ~ 5 x 10~^. k is 
the washout factor Q and e, = ±j:,^^{f{\^) + ^ith f{x) = 

V^M^)] and gix) = 2^. 

The recent experimental bound on r]B isj25[: 

VB = (6.5±°:^) X 10-1° (17) 

In the following sections, we will use these expressions to evaluate the baryon asymmetry, 
which along with present neutrino oscillation results will then give us a range of parameters 
in the seesaw formula. We will use these parameters to calculate the edms of the electron 
and the muon using the allowed range of SUSY breaking parameters mentioned above. 

We present three different cases. In at least one case, the effects are big enough that they 
are within reach of contemplated experiments for the elctronjl^]- Planned searches for the 
muon edm^J] will also throw important light on this issue. 



VI. FIRST MODEL FOR NEUTRINO MASSES AND ITS PREDICTIONS FOR 
LEPTON EDMS 

In this section, we consider the first of three examples for neutrino masses and obtain the 
predictions for lepton edms in the presence of the / terms. In this example, we take the 
and Yi to be diagonal and assume that all neutrino mixings to arise from arbitrary complex 
/ coupling^]. In the notation of sec. 2, this means that W is a identity matrix. Clearly 
therefore, in the absence of the / term contributions, edm becomes very tiny (< 10~^^ ecm. 
The / couplings lead to an enhancement of the effects of leptogenesis phases in the lepton 
edms as we see below. 

As for the Yiy, we consider two subcases: (a) in the first case we let the diagonal elements 
of be free and (b) in this second case, we let them be proportional to the charged lepton 
masses. 
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In both cases, we can write the matrix elements of / in terms of the diagonal elements of 
i.e. (ce^, e, l)y and the neutrino mass matrix. For this purpose, we take the light neutrino 
mass matrix to be: 

/ee" he"" de \ 



niQ 



(18) 



he"" 1 + ae 1 
\ de 1 1 + beJ 
Here (a, b, d, e, h) are order one coefficients, Am^^^ ~ Aitlq and the exponents {n, m) in 
the (1,1) and (1,2) entries have to be at least 1, but can be larger. This matrix provides a 
good fit to all neutrino data. Using these, we get for / matrix: 

/(a + 6)cV crfe^ -cde^ \ 



f 



V 



(19) 



/ 



V 



(20) 



cde^ —d^e^ dhe^ 
-cde^ dhe^ (e-/i2)eV 
We now calculate the edms and ris in this model. For example, at tan/3 = 30: 
we have / at 2 x 10^^ Q^y 

/-1.14xlO-^ -1.56x10-2 0.293 \ 

-1.56x10-2 0.379 - 0.353i -0.21 + 0.55i 
0.293 -0.21 + 0.55Z 0.058 - 0.091 * i / 

We now assume that the Dirac neutrino mass is proportional to the charged lepton mass 
and are 0.0337, 6.94, 123.398 GeV. Using these parameters, we find sin2 26'0 = 0.8, Am^ = 
5 X 10"^ eV^ and sm'^26A = 0.98, Am\ = 6 x 10"^ eV^. In the calculation we decouple the 
neutrinos at the respective mass scales. The baryon to photon ratio is t]b ~ 6 x 10^^°. We 
will use these inputs to evaluate the lepton edms. Our results for edms are given in Fig 1 
and 2. 

In Fig.l we plot the electric dipole moments of electron and the muon as functions of Aq 
and mi/2 and mo for tan (3 = 10. We do not show explicitly the values of mo since we choose 
their values in such a way so that the relic density constraint is satisfied in the only available 
stau-neutralino co- annihilation region for the parameter space. We apply the recent relic 
density constraint i.e. 0.094 < Qh"^ < 0.129 (2(T)j26l]. Using this constraint, the maximum 
value of mi/2 is found to be 800 GeV for tan/3 = 10. The magnitude c?^ can be as large as 
10-29 ecm for Aq = 800GeV, where as the de can be as large as 10-^^ ecm. 

In Fig. 2, we show the dfj_ and the de for the same model using tan/3 = 30. As expected 
the edms are larger in this case with dn as large as 4 x 10-^*^ ecm. 
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FIG. 1: Electric dipole moments of electron \de\ (solid line) and muon (dotted line) for different 
values of Aq for tan (3 = 10. 



VII. MODEL II FOR NEUTRINO MASSES AND PREDICTIONS FOR LEPTON 
EDM 

In this case, we choose Yi and / matrices to be diagonal and / equal to Diag{e^j,e^f, 1) 
but keep a general form for as follows: 

/ he^ de^\ 



Y„ 



he ae 

3 



be' 



\de 



be^ e'^ 



(21) 



For e ~ ej ^ 1 and a, 6 ~ 1 but complex, we get large solar and atmospheric neutrino 
mixings. We have done a detailed fit to the neutrino parameters keeping a, b, h, d complex 
and have calculated the baryon asymmetry and the lepton edms. For example, at tan (3 = 10: 
we have / at 2.5 x 10^^ GeV 



f = Diag[{6 Ax 10-^ 1.6 x 10"^ 1 )] 



(22) 
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FIG. 2: Electric dipole moments of electron de (solid line) and muon (dotted line) for different 
values of Aq for tan (3 = 30. 



The Dirac neutrino coupling is 

/ 1.67 X 10-2 - 1.43 X lO-^i -2.31 x lO'^ - 7.7 x lO'^i \ 

= 1.67 X 10-3 - 1.43 X 10-2^ 6.4 x 10^^ ^_ i i ^ iq-i^ 8.3 x 10"^ - 1.14 x 10~if23) 
V -2.31 X 10-2 - 7.7 X IQ-H 8.3 x 10"^ - 1.14 x IQ-H -0.57i J 

Using these parameters, we find sin^ 26^0 = 0.83, Ami, = 5 x 10"^ eV^ and sin^ 26a = 0.88, 
Am\ = 2 X 10~2 eV^. The baryon to photon ratio is ~ 6 x 10^^°. 

In Fig. 3 we plot the electric dipole moments of electron and the muon as functions of 
Aq and mi/2 and mo for tan/3 = 10. The relic density constraint is satisfied in the only 
available stau-neutralino co-annihilation region as before. The magnitude d^u can be as 
large as 10~28-5 ecm for Aq = 800GeV, where as the dg can be as large as lO-'^^ ecm. 

In Fig. 4, we increase the scale where the universality and proportionality assumptions 
are made to 10^^ GeV and study its effect on d^u as a function of mi/2 for Aq = 300 GeV. 
One needs to be careful about the size of the couplings at this new scale. The effects of the 
/ couplings in this new region increase the dipole moments. The top line shows the dipole 
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FIG. 3: Electric dipole moments of electron \de\ (solid line) and muon (dotted line) for different 
values of Aq for tan (3 = 10. 



moment for the increased scale and we conclude that the edm increases as the scale moves 
up. 

Note that all these values are considerably above the values in the absence of the / 
couplings [2^. Also note the violation of the mass scaling law for the dipole moments in the 
presence of the / couplings. 



VIII. THE CASE OF TWO vr'S AND THE 3x2 SEESAW 



In this section, we consider the predictions of a class of models with 3x2 seesaw considered 



in various papers ^, |27|, |28|. First point to note is that the general results we discussed for 



the case of three right handed neutrinos in sec. 2 do not apply to this case. In particular, as 
has been shown by Endoh et alQ], the leptogenesis phase and the PMNS phases are directly 
related in this model unlike the three RH neutrino case. 

To proceed with the investigation of this model, one can choose a basis such that and 
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FIG. 4: Electric dipole moments of muon dn for Mq = 2.3 x 10^^ GeV and 10^* GeV 



-lis 



/ are diagonal. The most general 3x2 matrix Yi, can now be parameterized as: 

/O m2 \ 
= V[ ]U 
VO mj 



(24) 



where V 



c s\ / e 



«7-R 







-s c/ V e-'T« 



and [/ = 023(^23)f/i3(^i3,<5)Oi2(ei2)PL where 







e-'^^ 



In order to calculate the edms and tjb, we will work with the horizontal model of Ref. [27! 
29|, which leads to a specific realization of the 3x2 seesaw formula. The model is based 
on the gauge group SU{3)c x SU{2)l x U{1)y x SU{2)h with fermion assignments given as 
follows: 



Table I 
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= (L„ L„) 


(1 2 -1 2) 




fl 2 -1 1) 


^ / c _ c\ 
1-j — VP ' ^ / 


fl 1 -2 21 




(1,1,-2, 1) 




(1,1,0,2) 


< 


(1,1,0,1) 


_ — V y 


(1, 1, 0, 2) 




(1,1,0,2) 




(1,2,+1,1) 




(1,2,-1,1) 


Ah 


(1,1,0,3) 



Table caption: We display the quantum number of the matter and Higgs superfields of 
our model. 

Here ive,/i,r denote the left handed lepton doublet superfields. Other symbols are self 
explanatory. We arrange the Higgs potential in such a way that the SU{2)h symmetry 
is broken by < Xi >=< X2 >= vm] < X2 >=< Xi >= and < A^.s >= v'h^ where 
VH,i,v'fj ^ Vwk- The vevs for x are chosen so as to cancel the D-terms and leave supersym- 
metry unbroken below the scale of horizontal symmetry breaking. 

The Yukawa superpotential for this model is given by: 

Wy = h,{L,H^vl + L^H^vD + /iiL.(z/^X2 + ^:Xi)HjM - z/iV^^rar ■ A^N' (25) 
^(LeX2 - L,xi)Hdr' + ^LrHa{ii'x2 + e'xi) + KLrH^r' + h'^{L,e'' + L^^^')Ha 

The parameters in the above equation have been determined to fit neutrino mass data at 
low energiesj2^. We use them in our calculation. We expect the predictions to be typical 
of most 3x2 models. From 0|, we see that at tan/3 = 30, the charged lepton mass matrix 
is given by (at the scale 10^'^ GeV): 

/ 7.5 X 10-3 -1.135\ 

Ml = 7.5 X 10-3 0.247 . (26) 

v2.14xl0-2 0.0984 0.85 / 
This determines all the Yukawa couplings responsible for charged lepton masses to be h'^ = 
1.135 X 10-\ h' = 1.29 X 10-3, h' = 1.466 x 10-^ and h' = 9.8 x 10-^ with K2 = 10 and 
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FIG. 5: Electric dipole moments of electron de (solid line) and muon (dotted line) for different 
values of Aq for tan /3 = 30. The model parameters are described in the text. 

^ = 4.588. We get the correct values of charged lepton masses at the weak scale from the 
above matrix. We use the MSSM RGEs between the horizontal scale and the weak scale. 
The / is given at the horizontal scale: 



Using the same ks as above, Hq = 9.2 x 10 ^ and hi = —2.92 x 10 ^, we find sin^ 26q = 0.87, 
Ami = ''•5 X 10"^ eV^, sin^ 29a = 0.88, Am\ = 2.48 x lO^^ eV^ and l^/esl = 0.11 which are 
within the experimentally allowed regions. We find 77^ ~ 6 x 10^^°. 

In Fig. 5 we plot the electric dipole moments of electron and the muon as functions of 



density constraint. We find that unlike other models, both dmu and dg are of same order 
and as large as 10~^^ ecm. 

To be completely phenomenologically consistent, we need to study the profile of lepton 
flavor violation for all these models. For model 1, this has been studied extensively in 3| 





(27) 



mi/2 and tuq for tan/? = 30 and 



— 



300 GeV. We again satisfy the relic dark matter 
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FIG. 6: The branching ratio for /x ^ e + 7 for model 2. 



and the results are within the present bounds but quite testable as has been emphasized 



Similarly, the branching ratio for /i — 67 in model 3 can be found in ref.j29|. We will 
therefore only calculate BR[yU —>■ 67] for model 2 for our choice of parameters and present 
them in fig. 6 for different values of mi/2 and Aq. We again demand the relic density constraint 
to be satisfied in this parameter space. We find that the most region of the parameter space 
is allowed by this branching ratio. The BR for r ^ /X7 is ~ 10^^ — 10~^^ in the same region 
parameter space. 



IX. SUMMARY AND CONCLUSION 



In summary, we have calculated the electric dipole moment of the electron and the muon 
in supersymmetric seesaw models where right handed neutrino masses arise from renor- 
malizable coupling in the superpotential. Using this we have attempted to probe the CP 
violating phases in the right handed neutrino mass matrix responsible for leptogenesis. We 
consider the minimal gauge group which allows for this. We make the simple and widely 
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used choice for the supersymmetry breaking parameters as in the mSUGRA models so that 
no new phases enter the discussion, other than the phases in the RH neutrino mass matrix 
which manifest at low energies and in leptogenesis. We find that due to the fact that the 
right handed neutrino masses arise via renormalizable couplings, the lepton edms get en- 
hanced and in some models come within the range accessible to proposed experiments. For 
instance, we can have de as large as 10~^^ ecm. The largest values for muon edm we predict 
are at the level of 10~^^-^ ecm. Also, the scahng law dmu/de — rrin/me does not hold in 
general. All the parameters of our model are chosen so that they are consistent with present 
neutrino data and the produce the required baryon asymmetry. The smallness of the edm 
values in most cases is due to the assumptions of complete universality of scalar masses in 
the MSSM and the proportionality of the A terms to the Yukawa coupling. To the extent 
that collider experiments have the potential to confirm or rule out the mSUGRA models in 
the near future, these results can teach us important things about the origin of matter as 
well as about the origin of neutrino masses. 

Clearly if the assumtions about SUSY breaking terms is relaxed, one can enhance the 
edm values further. One must then make sure that new phases not connected to leptogenesis 
are not responsible for the enhancement. In the minimal model there is no such confusion. 
This possibility is presently under consideration. 
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